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Introduction 

The  forced  oscillation  of  a  rigid  body  in  the  surface  of  a  fluid 
is  investigated  in  this  paper.  John  [1]  showed  that  the  problem  of 
steady  state  oscillations  of  a  body  in  a  free  surface  can  be  reduced  to 
a  Fredholm  integral  equation.  However  the  solution  of  the  integral  equa¬ 
tion  is  difficult  except  for  the  case  of  special  cross  section  geometry. 
Peters  and  Stoker  [2]  as  well  as  Haskind  [3]  developed  a  method  of  solu¬ 
tion  for  a  thin  ship. 

The  present  paper  concerns  the  three  dimensional  problem  of  a 
ship  with  small  draft.  VJe  consider  the  circular  and  elliptic  disks  and 
determine  the  dependence  of  the  added  mass,  added  moment  of  inertia, 
and  damping  factor  on  the  frequency  of  the  forced  oscillations.  A  body 
form  of  disk  provides  large  wave-making  effects  so  that  the  results  will 
serve  as  a  complement  to  the  thin  ship  theory.  Two  dimensional  aspects 
of  this  problem  have  been  treated  in  [4]  and  [5]. 

I,  General  Formulation 

We  suppose  the  half-plane  y  <  0  to  be  filled  with  an  incomprcss- 
ible,  nnviscid  fluid,  ^  0  corresponding  to  a  free  surface,  and  x 
and  z  denoting  rectangular  axes  on  that  surface.  It  is  assumed  that 
all  motions  of  the  fluid  are  irrotational,  time-periodic,  and  small 
enough  so  that  the  problem  can  be  linearized  by  neglecting  squared 
t  enns . 

V/e  suppose  that  the  fluid  motion  is  produced  by  a  ship  which 
is  placed  in  the  surface  of  the  fluid  at  rest  and  set  into  the  force- 
ed  periodic  oscillation.  When  transients  are  passed,  the  resulting 
fluid  motion  can  be  considered  to  bo  time-periodic  with  frequency 
The  irrotationality  implies  the  existence  of  a  velocity  potential 
0(x,y,z;t),  and  the  periodicity  in  time  means, 

(I.l)  0(x,y,z:t)  -  Re[U(x,y,z)  e'*^^]. 

Next  expressing  the  equation  of  ship  surface  in  the  equilibrium 


position  as. 


(1. 2) 


2 


y  -  S*(x,z), 

we  assiune  that  the  form  of  the  ship  S*(x,z)  satisfy  the  following  geo¬ 
metric  conditions: 

S*(x,z)  0  along  the  edge  of  the  equilibrium  water  plane 

area  which  is  bounded  by  the  curve  C, 

(1.3)  S*(0,0)  -  c,  maximum  draft  at  the  center, 

S*(x,z)  -  S»(-x,z),  and  S»(x,z)  -  S*(x,-z). 

In  view  of  the  symmetry  the  co-ordinates  of  the  mass  center  can  be 
written  as  (0,y^,0),  If  we  e^qsress  the  position  vector  of  the  mass 
center  as  ^  “  iX  +  3^  its  components  are, 

Ut)  -  Re(jr*  e'*®”^), 

(1.4)  7(t)-y  +  Re(7*  e*‘«^  ), 

Z(t)  -  Re(Z»  e'*«'^). 

V'e  call  7,  f,  and  Z,  the  surge,  the  heave,  and  the  sway,  respectively, 
and  X*,  Y*,  and  Z*,  the  amplitude  of  the  respective  motions. 

Now  let  us  introduce  a  set  of  moving  co-ordinates  (x' ,y* ,z'  ) 
whose  origin  is  attached  to  the  centroid  of  water  plane  of  the  ship  in 
the  equilibrium  position.  V/e  see  that  (x',y*,z')  will  coincide  with 
fixed  co-ordinates  (x,y,z)  if  the  ship  is  at  rest.  From  rigid  body 
dynamics,  the  velocity  of  a  particle  in  the  ship  at  any  time  is  given 

by, 

(1.5)  r*-R  +  ^x(r'  -Jy*) 

where  r  ••  ix  +3y  '♦'Rz,  r'”  '*■3?'  ,  <-0  =  iO  +JG  +kO  , 

X  y  z 


and 
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e^(t)  -  Re(G* 

(1.6)  e  (t)  -  Re(G*  e-‘'^  ), 

V  J 

G,(t)  -  Re(G*  ). 

z  z 

V;e  call  0„f6„,  and  G  ,  the  roll,  yaw,  and  pitcn,  respectively,  and 
X  y  z 

Gi,G*,  and  Q*,the  amplitude  of  the  respective  motions. 

X  y  z 

In  order  to  show  the  meaning  of  small  motions,  we  transform 
the  space  variables  by, 

(1.7)  X  -  kx,  y  -  ky,  z  -  kz, 

where  the  wave  number  k  is  k  - <<^/g  ■  2n/\  ,  \  being  the  wave 
length  of  free  waves  of  frequency The  normalized  amplitudes  for 
the  linear  motions  are  given  by  X*  =  oX’ ,  Y*  -  oY’ ,  and  Z*  -  oZ’  , 
o  being  a  small  parameter.  We  observe  that  a  small  means  the  ratio  of 
actual  amplitudes  to  the  wave  length  is  small.  Similarly  the  amplitudes 
for  rotational  motions  are  given  by  ©•  ■  oG^  ,  G*  -oGy  ,  and  0* 

Neglecting  squared  terms  means  that  the  potential  of  the  fluid 
motion  produced  by  the  forced  oscillation  will  have  the  form, 

(1.8)  0  =  a0', 

where  0'  and  its  derivatives  are  bounded  and  we  are  simply  neglecting 
terms  involved  in  all  our  formulations.  For  instance  we  express 
Bernoulli's  equation  as, 

(1.9)  P  -  -pgy  -  P0^(x,y,z:t)  =  -pgy  +  p  Re[i<iU(x,y,z)  e"*®’^], 
v/hile  the  surface  elevation  is, 

Y.(x,z:t)a;-  i  0^(x,O,z:t). 


(T.iO) 


k 


¥e  shall  now  consider  tho  boundary  conditions  governing  the 
potential  function  U(x,y,z)  in  (I.l).  As  the  fluid  is  assumed  to  be  in¬ 
compressible  and  Irrotational, 

-  0, 

or 

(1. 11)  ■0  in  y  <  0,  outside  the  ship. 

Neglecting  terms  in  the  free  surface  condition  is, 

^tt  * 
or 

(1.12)  U-  -  kll  -  0  on  y  ■■  0,  outside  the  ship. 

As  no  flow  occurs  across  the  surface  of  the  ship,  the  kinematic  condition 
can  be  written  as. 


0^  •  n(x,y,*:t)»£(x,y,z:t)  on  the  iimersed  surface 

y  -  S(x,i:t) 


where  £  is  the  normal  vector  to  the  surface  given  by, 

II  ■  i  cosCn^x)  +  J*  cos(n,y)  +  k  cos(ri,z). 


and  £  is  the  velocity  vector  of  a  point  (x,y,z)  on  tho  surface.  We  re¬ 
mark  here  that  the  immersed  surface  in  motion  S(x,z:t)  differs  from  the 
equlibrium  surface  S*(x,z)  introduced  earlier.  However  it  is  consistent 
with  the  omission  of  terms  in  to  require  this  condition  to  be  satis¬ 
fied  on  the  equilibrium  surface. 

If  we  denote  tho  velocity  potential  for  surge,  heave,  sway,  roll,  yaw, 
and  pitch  by  0^  J  =■  1,2,... 6,  respectively,  the  linearity  of  the  pro¬ 
blem  permits  to  write  the  total  potential  as  0  “  0j  •  Furthermore 
the  components  satisfy, 


-  Aco-.i  n,  c). 


Yror,ln,y), 


(0^)n  '  Zcos(n,z), 
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(0^)n  ■  0jj[(y-yo)cos(n,z)-zco3(n,y)], 

(0j)n  ■  0y[zcon(n,x)“xcos(n,2)], 

hence, 

— i<l*co3(n,x), 

( — i<^e^C  (  y-y*  )co8(n ,  z)  -zco8(  n  ,y)  ] , 

(1.13) 

(Uj)^  — i<<9yCzco8(n,x)-xco8(ri,2)] , 

(U6)jj  — i^0j[xco8(n,y)-(y-y,)co8(n,x)],  on  y  -  S(x,zit). 

Finally  at  large  dlatance,  the  propagating  dlaturbanee  muat  have  the  form 
of  a  radially  outgoing  progreaalve  wave,  that  la, 

(1.14)  U(x,y,z)  •  O(^)  as  r  4e>, 

where  r  -J  x*  +  z*  ,  and  0  -  arctan  (4)  . 

To  show  clearly  the  dependence  of  the  solution  on  parameters,  we 
now  introduce  the  dynamic  pressure  functions  Pj(x,y,z)  J  -  1,2,... 6, 

by. 

g  ^  Pj^(x,y,2)  -  i^Uj^  (  p  ^  ), 

R  ^  P2(^»y»*)  "  ^  k*  f 

g  f  P3(x,y,z)  “  ^*^^3  (  f»  k»  f 

(1.15) 

g  O*  ^  p^^(x,y,z)  »  P  7  )» 

g  O*  ^  p^(x,y,z)  -  Uj(  ^  )» 


t  f*  k»f 


then  the  conditions  (l.ll)  -  (1.14)  can  be  expressed  as. 
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(1.16) 

(1.17) 


a.,  _  0 


V  P 


in  y  <  0,  outside  the  ship. 


(1.18) 


(Pj)y  -  (Pj)  *  0  on  y  •  0,  outside  the  shipj 

(Pl)n  "  co3(n,x),  (P2)n  "  co®(ll»y)*  “  cos(n^z), 

(p^)n  “(y“yo)co8(jl>2)’2co8(n,y), 

^^5^n  "  zcos(n,x)-xcos(n,z),  (p^),,  “Xcos(n,y)-(y-y,)cos(n,x) 


(1.19)  Pj  ^  e  -  0(i) 


on  y  ■■  S(x,z:t)^ 

as  r  4  CO. 


Equations  (1.18)  can  be  rewritten  as, 

-S 


<Pl^n 


^.52 

'  X  z 


,  ^  -(y-ye)s  -2 

(p.) - ^ 


^  "  /l  +s*  +s* 


2  ' 

X  z 


(1.20) 


JT+S*  +S^  * 

'  X  z 


•zS  +  xS 


2  +s*  * 


(Pq).,  - 


iPA)«  - 


x  z 

x*-(y^c)s 


^  +5^  +S^  *  ^  "  yi~+S*  +s* 

»  X  Z  y  yc  r. 


X  z 


on  y  -  S(x,z:t). 

Here  if  we  neglect  terms  in  o^,  the  quantities  S(x,z:t)  in  (1.20)  can  be 
replaced  by  S*(x,z). 

By  (1.15)  we  have  Bernoulli's  equation  in  the  form. 


?!  ’’-Pr.  k  f  Rc[pj^(x,y,z)  e**^]. 
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(1.21) 


Pg  — Pf5  f  PK  1“  Re[p2(x,y,z)  e'*'*], 

P3  —PR  k  f  R«CP3(*»y»®)  ]» 

P^  —PR  f  ♦  P^*  f  Re[pjx,y,z)  ], 

P5  —PR  f  +  P^  f  Re[pj(x,y,z)  e’*''*  ], 


II.  Shallow  Draft  Approximation 

In  this  section  we  relate  the  first  order  forces  and  moments  ex¬ 
erted  on  a  ship  which  oscillates  with  six  degrees  of  freedom  to  the  added 
mass,  added  moment  of  inertia,  and  damping  factor.  Then  wc  develop  the 
perturbation  procedure  for  a  ship  of  small  draft. 

The  three  components  of  the  fluid  force  relative  to  the  fixed  co¬ 
ordinate  system  are, 

••//gP(x,y,z:t)co8(n,x)  dS, 

(II. 1)  Py  -//3P(x,y,z:t)co8(n,y)  dS, 

Fj,  - Jj|^gP(x,y,zst)cos(n,z)  dS, 

where  P  is  the  pressure,  and  S  here  represents  the  immersed  surface  in 
motion.  The  three  components  of  moment  relative  to  the  fixed  co-ordinate 
axes  are, 

-•/jgP[(y-7)cos(n,z)-(z-Z)cos(n,y)]  dS, 

Gy  gP[(z-Z)co5(£,x)-(x-X)co3(n^,z)]  dS, 

G^  =y^^P[(x-7)cor!(£,y)-(y-7)coD(£,x)]  dS, 


(11.2) 


Prom  the  calculated  results  of  forces  and  moments  presented  in  [1],  for 
the  ship  satisfying  the  s^naetry  conditions  given  by  (1.3)  we  find, 

(II. 3)  Fy  dS  +  Kg  -  opgf’A  +  0(o*), 

%  "■“‘’//s'^t  ^  *  0(a*), 


and 


°x  ■*p//8*^t^®®°®^E^y)“^y“ye)co®(jl*®))dS  -  apg0^//^(*‘'4^)dxd»  +  0(o*) 


(II.4)  Oy  ^pff^,^\lxcoB(ntz)-ieoa(n,x)'i<3S  +  O(o^), 

"«pZZ8«<^t^(y“ye)«o*(ll>x)**c‘>®(lI*y))‘*S  -  opg©^//^(x*4j5drd£  +  0(a*) 

i^ere  S*  denotes  the  immersed  surface  in  the  equilibrium  position,  and 
A,  water  plane  area  of  the  ship. 

Next  we  express  the  forces  and  moments  for  the  oscillatory  motion, 
up  to  teims  of  order  a,  in  the  dimensionless  form, 

r  Re(Fxe'*'^  )»  -Re[  jf^,p^coa(T^x)dS  e"‘*^  ), 

~^gY'  MFyO'**^^)-  -Re[yZg.P2C08(n,y)dS  e'‘**  ]+  J^T«e(e*‘*^  ) 

-k*ARe(e-‘^^), 

^^^1  Ro(F^e*‘^  )•  Hlo[ yZg.P3Cos(n,z)dS  ], 

and  ^ 

—jjgr  Re(0^e‘*'^  )-  Re[  y^,p^(zcos(n,y)-(y-y,)co3(n,z))dS  e'‘^^  ] 

”  Re(e'‘*’^^), 

(II. 6)  RR(rfyC'‘*^^)-  Re[yyg,Pj(xcos(n,z)-zcos(in,x)d:j  e***^^], 
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apgae, 


^  Re(0^e‘*^*)-  Re[  //g.p^((y-y,)co8(ji,x)-xcos(ii,y))<lS  e**^  ] 

-  +  §  )dxd8  Re(e**^  ). 


Now  we  seek  to  express  the  oscillatory  forces  and  mcments  as, 

Re(V‘‘^^)“  ■^x^'^x^*  Re(G^e-‘'^)- 

(II.7)  Re(Fye**'^)-  -My?-iIyf^flgRe(e"**^  )-pgK,  Re(Gye‘‘'^  )-  -lyS^y-flySdy, 

Re(Fj^e*‘«*  )-  Re(G^e  )-  -1^5^ 

where  K  and  1  are  called  the  added  mass  and  added  moment  of  Inertia,  re* 
sepectively,  while  fl  and  il  are  called  the  damping  factors  for  translation 
and  rotation,  respectively. 

The  substitution  of  (l*/»)  and  (1.6)  in  (II. 7)  yields, 

Re(  F^e  )-  oReC  (d*M^+ldN^)|'  e  ***< ) , 

Re(F  e*‘'^)-  aRe[(<%+idN^)J*e***^  )+MgRe(e*‘'^0-apgRe(?’e*‘'‘  )A, 
y  /  /  K 

Re(F  e**^*)-  oRe[(d^K  +idN  )f*e  *»»*], 

Z  Z  Z  K 

(II.0) 

Re(G^e’‘®*)-oRe[(d^I^+i^^)^  ©j^e'***  ]-ap|^e(ej^  ^  e’***  )  //;^(**+|‘  )dx'lz. 


Re(G  e-i^*)-aRe[(4''l  e'e'*^^  ], 

y  j  j  ^  j 


Re(G^e*i^*)-oRe[(<J^I^+idH^)^  e^e**»*  l-<ipRRe(©i  f  e-‘«*  )  ff ^{x^^)dx6l. 


hence  in  the  dinenni onlesn  form. 


^  i!0(  v“'“‘  '■ 
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^  R.<Py.-«'  )-  k«R.t(  +li)  1*^  H.(,-«'  ) 

-k*ARe(e-‘<^^  ), 

^  ReCF^e-**^)-  k»R«[(  ], 


I  ff 


(II.9)  Re(G^e*‘*^  )-k^e[(  +i^)  e-‘«*  ]-  i  //^(**-»|’)dxd*Re(e-‘»<  ), 

u4  I 

Re(0  )-k<Re[(  ^  +i-^)  «•*«*  ], 


otpga^  '  7 


P  Ptf' 


J5i^R.(0,,-‘«)H,''Rrf(  .-i«]-i//j(:t“*f)dxddl.< .-•*') 

& 

Therefore,  equating  (II. 5)  and  (II. 6)  to  (II. 9)  we  obtain, 

k*ffyp  -//^  ,(pj^)j.co8(n,x)dS,  k*  Rj/Pd— //e  •  ^  Pp  i  "  *^Ii**^ * 

k*My/p  — ^.(pppC08(n,y)dS,  k*Hy/pd— //^,(pp^co3(n,y)dS, 

k*M^/p  ff^,{p\)j,co9injz)dS,  k'flypd— //g,(p^)jC08(n,*)dS, 

(II. 10)  k^r^p  -/^.(pJ^)pCBco8(n,y)-(y-y,)co8(n,z))dS, 

k^Hjj/pd-  y/g,  ( P^  if  *co8(  n,y )  -( y-y*  )coa(  n ,  z)  ]dS , 

Vp  -//s  ,(pj)^Cxco8(n,z)-zco3(n,x)]clS, 

k^Hy/p4-  //g, ( Pj ) J [ XC03( n , z ) -ZC08( n , X ) ]  dS , 

k^f^/p  -^,(p^)j.[(y-y*)cos(n,x)-xcos(n,y)]dS, 

k^H^/p<S"  ^.  ( p^  t  ( yy*  » 
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where  the  pressure  function  p(x,y,*)  is  resolved  to  the  real  and  imaginary 
fart.  a.  Pj  -  (Pj)^  *  l(Pj)i  j  . 

We  proceed  to  consider  the  shallow  draft  approidmation.  First  we  re¬ 
mark  that  the  smallness  of  draft  means  the  ratio  of  actual  draft  to  the 
wave  length  is  small.  Therefore,  all  quantities  are  to  be  developed  as  a 
power  series  in  e  ■2Ke/K . 

Let  the  profile  of  a  ship  of  small  draft  be, 

(11. 11)  y  -  eS'(x,z)  “a  <  X  <  -b  <z  <  b. 

Here  s’(x,z)  satisfies  the  following  geometric  conditions! 

S*(0,0)  -  1,  S’(x,z)  -  s'(-x,z),  S’(x,z)  -  S’(x,-z). 

Since  -  1  -  f^C(s;)Ms;)‘') 

nip2 

we  obtain, 

5^ .  (..s; .  o(.“))  52  *  [  1  «(.>)]  jf  M-es; .  jf . 

If  we  call  the  right  hand  members  of  (1.20)  Rj(x,z)  J  -  1,2,... 6, 
these  can  be  expanded  in  the  form, 

R4(x,z)  ==•  c"  r"(x,z). 

The  first  few  terms  are, 

H*(x,z)  "  0,  R^(x,z)  -  1,  R^(x,z)  -  0,  RJ^(x,z)  — z, 

(11.12) 

R*(x,z)  -  0,  R2(x,z)  -  x,  ...  . 

Next  let  us  assume  that  the  pressure  can  be  expanded  in  the  form, 

P,(x,y,z)  -  2Z  ‘="'p'!(x,y,z) 
mrO 


(II. 13) 


J  -  1,2, ...6 
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By  Taylor's  theorem  the  function  Pj  can  be  expressed  as, 

p'![[x,eS'(x,*),z]  -21^  teS’(x,z)]"  p'J(x,0,8). 

J  n-0  dy"  J 

Now  we  can  write  the  left  hand  members  of  (1.20)  up  to  order  e  as, 

(II. H)  [  1  +  0(e*)]  2Z^  e"**‘"[s’(x,*)]"  Pj(x,0,s) 


n' 

♦  C-eS'  +  0(c*)]  L  -m+n 

*  i-6 


n+l 


dwy 


-  RJ(x,s)  eR^(x,B) 


on  y  -  cS’(x,b). 


Equating  the  coefficient  of  the  like  powers  of  e  we  obtain, 

dp*(x,0,») 

dpUx,0,a)  dp!(x,0,z)  dp*x(x,0,z)  dV.(x,0,s) 

— J - -  R’(x,8)  +  s’  — - +  s'  - s' - J -  . 

dy  ^  dx  *  dz  dy* 

Kore  generally  the  boundary  condition  for  Pj  is  of  the  form, 

dp'?(x,0,s) 

'  ■  ■  Rj(x,b)  +  within  the  water  plane  of  the  ship, 

n  k 

Here  Ej  consists  of  Pj  for  k<  n>l  and  their  derivatives  evaluated  at  y  ■ 


Therefore  Pj  can  be  determined  recursively  by  solving  boundary  value  pro¬ 
blems  of  the  following  form; 

Find  a  function  w(x,y,z)  such  that. 


(11.15) 


-  0 


in  y  <  0,  outsidr  the  ship. 


(11.16) 


w  -  w  -  0 

y 


on  y  ■  0,  outside  the  ship. 
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(11.17)  Wy  -  h(x,z)  on  y  -  0,  within  the  hull  of  the  ship, 

(11.18)  w  -  e**^  e^**  -  0(i)  as  r  ♦oo. 


Note  that  h(x,z)  is  a  given  function,  and  for  w  ■  p*  in  particular  it 
is  given  by  (11.12). 

V/e  again  make  a  change  of  independent  and  dependent  variables  by 
introducing, 

x-x/a,  y-y/a,  ?-z/a. 


(II *19)  w(x,y,z)  ■  a  V(x/a,  y/a,  z/a)  where  a>ka  . 


Then  the  boundary  value  problan  becomes: 
Find  a  function  V(St,y,z)  such  that. 


(11.20) 

•  0 

in  y  <  0, 

outside  the  ship. 

(11.21) 

Vy  -  a  V  -  0 

on  -  0, 

outside  the  ship. 

(11.22) 

Vy  -  h(S,z) 

on  y  "  0» 

within  the  hull  of  the  ship. 

(11.23) 

V  -  ^  c*? 

-  O(^) 

as  r  ^-oo  . 

Hero  a  coimnent  on  the  change  of  variables  is  in  order.  Since  x  kx  and 
a  ■  ka,  the  x  corresponds  to  the  ratio  of  actual  co-ordinate  x  to  the  half 
length  of  the  ship  a  .  Thus  in  the  new  boundary  value  problem  all  length 
dimer  Fslons  have  been  made  dimensionless  by  dividing  with  the  half  length 
of  the  ship  instead  of  the  wave  length. 

Now  we  turn  to  the  determination  of  the  added  mass  or  added  moment 
of  inertia  and  the  damping  factor  for  a  ship  of  small  draft.  Since  we  now 
have  y,-  0(c),  cos(n,x)  -  -cS^  +  O(c^),  cos(n,y)  ■  1  +  O(c^),  and 
cos(n,z)  -  -cS^  +  0(c‘-’),  the  following  expressions  can  be  found  from  (II. 10). 


u 


Note  that  for  convenience  we  adopt  the  notation  Pj(x,y,z)  for  previously 

introduced  pl(x,y,z).  S*  is  In  this  case  a  rerion  in  the  x-z  plane, 

J 

k*  M^p  —  //g.  [P3^(x,y,z)]^co8(n,x)dS 

♦e[p^(x,z)]p\eSj^dS-  0(e), 

similarly  — /f3.[p3^(x,y,z))^co8(n,x)dS-  0(e), 

k*#iy/p  — ^g,CP2(x*y»*)3pC08(n,y)dS 

"*  ♦e[pJ(x,z)]p}dS--  //^.[p*(x,z)]^dS. 

similarly  k*<Iy/ptf  — ^g.[p2(x,y,z)]j^co8(n,y)dS— ^,[p5(x,z)]  ^dS, 

k*fl^/p  —  //g,CP3(x,y,z)]^co8(n,z)dS 

- //a.{[p5(x,z))^+eS«^^2ijJ2i:  +e[p>(x,t)) JcSjdS-  0(e), 
slmiUrly  k*ll^/p<5  — ^.CP3(x,y,z)]iC08(n,z)dS-  0(e), 


(II.2i») 

k*I^p  -  ^^,[p^(x,y,z)]^[zco8(n,y)-ycos(n,z))dS 

-)5(,.[p];(x,z))r(  z+eS’eS;)dS.^^.  z[p*(x,z)]^dS, 
«^g,[p^(x,y,z)]j[zco8(n,y)-yco8(n,z)]dS»  z[p*(x,z)]j^dS, 

k^i  y/p  -  y7g,[p3(x,y,z)]^[xcos(n^z)-zcos(n^x)]dS 

“ //g.[P5(x,z)lp(-xcS’+zcS')dS-  0(e), 

k^H y/p^  - ^^,[p^(x,y,z)]^[xco8(n,z)-zco8(n,x)]dS-0(c), 


k*i^/p  -  //3,[p5(x,y,*)]j.tyco8(n,x)-xco8(n,y)]dS 
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-//a.Cp2(x,e)lr(-eS‘eS^-  x)dS— /^,  x[pj(x,z)]^dS, 

-//a,Cp^(x,y,z)]j[yco8(n,x)-xco8(n,y)l(iS— /j(^,x[pg(x,*)]^dS. 

Hence,  In  the  dlmenelonlees  form  we  obtain  the  added  maas  and  damping 
factor  for  heave  aa. 

My-  My/M»  .  -  i//5tp|lSr.S)l,MS, 

(11.25) 

Ny.  fiy/p9<r.  -  j//,(p5(Sf,5))i<is, 
the  added  moment  of  inertia  and  damping  factor  for  roll  as, 

V  “  a/4  ®  tp*(x,z)]j.dS, 

(11.26) 

ffypaV-  ?  [p|^(x,z)]j^dS, 

and  the  added  moment  of  inertia  and  damping  factor  for  pitch  as, 

I*-  ■  a//?  * 

(11.27) 

»*-  H/pI^5-  -x/fsr*  [pj(x,z)]jdS, 

where  is  the  im.ige  of  S*  under  th  '  transform.ition  7  x/a  ,  z  -  z/&, 

Ry  (11.19)  we  write,  Pj(x,y,z)  »  aV'^(x/a,y/a,z/a)  J  -  2,4,6, 
then  from  (11.22)  we  find, 

(11.28)  Vi(x,0,z)  «•  1,  vl(x,0,z)-  z,  v5(x,0,z)-  x  on  S, 

y  y  y 

nnd  from  (11.21)  we  h;ive. 


(11.29) 
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vJ(5r,0,2)  -avJ(x,0,?)  -  0  J  -  2,/*,6,  outside 


Ill .  Integral  Representation. 

We  will  discuss  here  the  integral  representation  of  the  solution 
of  boundary  value  problems  for  a  surface  obstacle  of  negligible  draft. 
Consider  a  region  bounded  by  the  surface  of  the  ship  and  the  free  sur> 
face.  We  shall  show  that  for  some  function  f(>r«2)  defined  over  the 
surface  of  the  shipi^the  potential  at  any  point  (x,y,2)  in  the  region 
is  given  by, 

(111.1)  ns,y.t)  -  ntV  ilo?  • 

Here  represent  the  Green's  function  given  in  [1]  eva¬ 

luated  at  0,  which  can  be  expressed  as, 

(111.2)  0(i!,jF,r,|,0,tl)  ■  V'* 

/CO 

where  the  Integral  sign  is  to  be  understood  as  integration  along 
the  positive  real  axis  except  for  an  arc  in  the  lower  half  plane  to 
avdid  the  positive  real  root  p  >  a  of  the  denominator.  Also, 

7  ■  •/  ,  and  ^  »  hence  R  denotes  the  distance  from 

a  point  (x,3P,?)  in  the  regio.^.  to  a  point  (|^,0,^)  on  the  surface  of  the 
ship.  Jg(Pr)  is  the  zero  order  Bessel  function  of  tho  first  kind. 

From  (III. 2)  we  see  that, 

(111. 3)  C--aG-i||u  in  jF  <  0. 

Now  we  define  tho  following  integrals. 


W(x,y,7)  -  —  //-  Hl,l)  , 
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L(5?,y,z)  -  i^jj^  H(3J,y,z)  d^d?:, 

where  ^ 

H(x,^.z)  -  J^Or)  e^y  ^ 

It  follows  from  (III.l)  that  V(x,y,z)  -  W(x,y»z)  +  L(x,y,z). 

For  a  function  f(xjz)  continuous  on  S  the  intef;ral  L  satisfies  the  condl* 
tion  (11.20)4  while  W  L  satisfies  the  condition  (11.23). 

From  (III. 3)  we  have^ 


V^(x,y,z)  -  aV(Sr,y,z)  -  VJ-(x,y>z)  in  y  <  0. 

^  w 


By  a  theorem  of  the  potential  theory, 


lim 

f  *0  ^ 
y  <  0 


nx,z) 

0 


on  S, 
outside  S, 


hence, 

(III./*)  V-(3r,0,z)  -  aV(x,0,z)  -  f(x,z)  on  S, 

(III. 5)  V-(Sr,0,z)  -  aV(x,0,z)  -  0  outside  S. 

y 

From  (11.22)  it  can  be  seen  that  the  potential  V(x,y,z)  given  by  (III.l) 
will  be  the  solution  of  our  bowidary  value  problem  if  f(x,'2)  is  chosen 
as  a  solution  of  the  following  integral  equation. 


h(}c,z)  -  aV(x,0,z)  -  f(X,z), 

or, 

(III.6)  f(x,z)  +  C(x,z,?,c)  d^dc-  h(x,z). 

Observe  that  the  zero  draft  approximation  leads  to  the  kernel  being  0 
itself  rather  than  the  normal  derivative  of  C  as  in  [1]. 

Next  we  calculate  the  kernel  of  the  integral  equation  explicitly. 


From  (III. 2)  we  have. 
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eoJ^(8F) 


(III.7)  G(Sr,0,a,|,0,?;)  -  I  +  2a\(r  -2-. 

Jo  B-i 


y  +  2a 


r<o  /■ooJjpr) 

\l/--c —  J  (pr)  dP  +  2a*\(/-  -2 - dp, 

^tp^*a*  X  p*“a* 


Let  us  express  (II I. 7)  as. 


(III.8)  G(3?,0,z,|,0,15)  -  I  4-  2alj^  +  2a2l2, 


where 


/«  PJ  (P?) 

I  ■  "■  ~  -  dp,  and  I^- 

^  /o  (p-a)(p+a)  ^ 


fOi  J  (p?) 

-dr - 2 - dp, 

L  (p-a)(p+a) 


’•'o  have 


r®  pJ  (p?)  ra-c  pJ  (pr;  f 

(III. 9)  •vlF  . .  dp  ■(  . . ■■■■—  dp  +  iyJ A ar)  *  f 

/o(p-a)(p+a)  /o  (p-a)(P+a)  Ja 


a-c  pJ  (p?)  ,  /•  oo  pJ  (p?) 


,a-e  pH"XpP)  r  co  pHjXp?) 

-  iyIJa?)  +  Re[  - - - dP  +  I  - - - dp], 

/  o(p»a}(p-fa)  /a+e(p-a](p+a) 
for  small  positive  e.  _ 

/CO 

Next  let  us  consider  the  integral  pF(*’XpP)/(p-a)(p+a)dp  inte- 

'  o 

grated  along  the  real  axis  except  for  an  arc  nwing  above  the  root  p  ■  a, 
then  we  obtain. 


'a+e(p-a)(p+a) 


(III. 10) 


-®  pH®(p?)  a-c  pH“‘(pr) 


(P-a)(P+a) 


(P" 


i)(P+a) 


►lSH«;(.r)  *  / 

'ft* 


phIXp?) 


'a+€(p-a)(p+a) 


a-c  pll‘’\pr) 


o  pH*|\pF) 


o(p-a)(p+a) 


a+c(p-) 


i)(P+a) 


From  (I II. 9)  and  (I II. 10)  we  find. 
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fCO 


(III. 11)  -  i|  J^(aPO  +  Rc[i|  Hj*(ar)  +A 


)o(p-a)(p+a) 


dp] 


■  ij  ^^(ar)  “  ^  +  Re  ^ 


00  pH^Vp?) 


dp. 


(p-a)(p+a) 

Here  we  can  deform  the  path  of  integration  into  the  positive  imaginary 
axis  by  setting  p  -  i^,  then  Ij^  becomes, 

roo  ’•/Hjjd^r) 

t 


-  _  /CO  7tr'ViYi 

I,  -  i|  J^(a?)  -  I  Y  (ar)  +  Re  --P---- 
1  2  o  2  o  jo  ^  ^  t 


-  4  J  (ar)  -  i  T  (ar)  -  Ratf  /  - i-S^  dY  1. 

2  0  2  o  '  i  ‘ 

Observing  that  the  integrand  is  purely  imaginary  we  now  express  Ij  as, 

(ni.i2) 


II"  |C  iJQ(a?)  -  Y^(ar)  ]. 


If  the  path  of  integration  runs  below  the  root  p  ■  a,  for  I2  we 


obtain. 


a-c  H“j(pr)  r  ®  H^Cpr) 


/a-c  n  ipr/  r 

(III. 13)  Ip  -  i^  J  (ar)  +  ReCj  - 2 - ap  +  / 

(p-a)lp+a)  'a 


a+e(p-a)(p+a) 


dp]. 


Repeating  the  same  process  employed  to  obtain  (III.IO)  we  find. 


_  _  H”\pr) 

(ni.U)  I,  -  Ijj  VaR)  - 


V/hen  the  path  is  deformed  into  the  positive  im.Hglnary  cixis,  I^  becomes. 


I  “  1- 

2 


«  n  P  r“ 

r-  J  (ar)  -  Y  (ar)  -  Ke[—  I  -  dV  ] 

2a  o  '  2a  o  n  ^ 


Therefore  by  a  result  in  L7]  we  find. 
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(in. 15) 


2a  o 


2a 


.^[iJ>?)  -S^(a?)]. 


Finally  from  (III. 12)  and  (III. 15)  we  obtain, 


(III.16)  G(2,2,|,1|)  -  I  +  2a  I  [iJ^(a?)-Y^(ar)]  +  2a*5|  [lJ^(a?)-S^(ar)] 

-  »  -  na[Y  (a?)  +  S  (ar)  -  i2J  (ar)], 
r  O  O  Or 


where  Y  (a?)  denotes  the  aero  order  Bessel  function  of  the  second  kind, 
and  S^(a?),  the  aero  order  Struve  function,  respectively.  Thus,  the 
kernel  of  the  Integral  equation  can  be  evaluated  explicitly  and  indeed 
this  is  another  reason  for  introducing  the  shallow  draft  approximation. 
Here  we  observe  that  as  ar  S^(ar)asY^(ar),  hence  the  Green's 

function  given  by  (III. 16)  becomes. 


G(x,a,](,?)a  I  -  2RatY^(a?)  -  J^(ar)] 


2  „  „r  8in(a?-itA) 


2  .  j  2na'  i(a?-n/U) 

7  ^  V ~7~  ® 

According  to  the  Fredholm  theory  the  integral  equation  (III. 6) 
will  be  soluble  if  the  corresponding  homogeneous  equation,  that  is, 

(III. 17)  f(x,2)  +  ~  11^  f(?,C)C(?,54,C)  d"?  -  0 

has  only  the  trivial  solution,  f*(x,a)  ■  0. 

It  has  been  shown  in  [1]  that  if  f*  ia  a  «^olution  of  (111.17)  then 

V(x,y,2)  -  f(LC)  G(x,y,a,^,0,C)  dC  , 

vanishes  identially  in  y  <  0. 
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Since  V*(5t,y,K)  vanishes,  by  (III, 4)  f*(x,z)  must  also  vanish.  Hence  the 
integral  equation  (III. 6)  is  soluble,  and  we  shall  have  a  solution  for 
our  boundary  value  problem  in  the  form  of  (III.I). 

When  the  solution  of  the  integral  equation  f(3[,z)  is  found,  we  want  to 
determine  the  added  mass  and  added  moment  of  Inertia  as  well  as  the  damp- 


Ing  factors.  By  substituting  (11.28)  in  (III. 4)  we 
related  to  the  density  function  f(x,z)  as. 

find  the  pressures 

(111.31) 

p*  -  aV^(x,0,z)  -  1  -  f^(x,z) 

for  heave. 

(III.32) 

p*  -  aV^(x,0,z)  -  z  -  f^(x,z) 

for  roll. 

(111.33) 

p^  -  aV^(x,0,z)  -  X  -  f^(x,z) 

for  pitch. 

where  f’^(x,z)  -  f^(x,z)  +  i  f^(x,z)  J  -  2,/,, 6, 


Therefore  we  obtain  from  (I II. 31)  and  (11.25), 

(111.34) 

l/Zs  ^i<5c,z)  dS 

from  (111,32)  and  (11.26), 

^x=  a//s  z-f5:(x,z)]  dS, 

(111.35) 

”x  “  •  a//s  * 

and  from  (III. 33)  and  (11.27), 

^7.  -  -l//s  X  -  f^(x,z)l  dS, 

«z"  l//s 


for  heave. 


for  roll. 


(111.36) 


for  pitch 
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IV,  Devftlopment  of  Asymptotic  Solution 

The  inapection  of  Green's  function  Riven  by  (III,]6)  shows  that 
it  is  a  function  of  a  parameter  a,  hence  the  solution  of  the  integral 
equation  (11 1. 6)  must  depend  upon  a.  Vie  develop  asymptotic  solution  for 
small  a  in  this  section  prior  to  treating  a  numerical  procedure. 

The  asymptotic  development  for  the  solution  of  two  dimensional  problem 
is  presented  in  [6], 

The  functions  appearing  in  (111,16)  have  expansions  of  the  follow¬ 
ing  form. 


oi  /  \m  M  2m  00 

J„(ar)  -  21  ^  %(?)  »  » 


2  ^ - (ar)^"’^-  -  ^  B  (?) 

"  m-O  [l,3,,.(2m-H)]*  ”  m-O  "* 


(IV,1) 


Y^(ar)  f  (log  ^  +  !r)  J^(a?)^  ^  (-l)"’^"(  £.  J 

m«l  n“l 


(ml)' 


—  2m 

(^) 


«  ^  [log  a  >'  A  (r)  * 

m-0  m=*0 


Hence  (1II,16)  can  be  expressed  as. 


(IV, 2)  G(x,z,I',^:a)  -  ^  “  2[log  ^  ^ 


00 


2m+l 


m=0 


m 


m»0 


m 


where, 


2m+2  .  2ni+li 

+  ^  13  (r)  a  -  171  2\  A„(r)  a  ], 

m=0  m=0 


^1  °  “  2*  ^2  "  *  ***  ' 


r^ 


13  -r,  B.  “  -  ,  -  ,  13 


r5 


1,3  *  2  1,3,5 


-  S*  log  ^  ,  Cj  ■=  -  (3f  +log  ^  -1  ) 


Rearran^ng  power  series  we  may  write. 
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(IV.3)  aG(x,2,|,?;;a)  -  2  «*  +2  21  log  a 

n-0  "  n-0  " 

with  new  coefficients, 

X  r 

=»  ^  ,  Oj^  "  in  “  iT  “  log  a2“”r,  •••> 

r^ 

Pq  "  “1»  Pj^  "0*  ^2  ^  *  ^3  ***  * 

Now  we  proceed  to  develop  the  asymptotic  solution  of  the  integral 
equation  in  the  form, 

00  CO  .  . 

(IV. 4)  f(x,2)  CO  2Z]EZ  loga)-^  . 

i-0  J-0 

Note  that  the  power  product  a^(a  loga)"^  can  be  always  ordered  as  to  their 
rate  of  vanishing  as  a  tends  to  2ero,  that  is, 

11m  *  IPg^A  .0  if  i»+J'  >  i  +  J,  or  i'+J*-  i  +  J  and  j>J' , 
a^(a  loga)^ 

Substituting  (IV, 3)  and  (IV. 4)  in  the  integral  term  of  (ill. 6)  we  obtain, 
//s  G(5:,*,?,^:a)  d^d^^ 


"  ’^11 B 


00  GO 


00 


i»0  J*0  n"0 


[  f  J  ( loga  )  ■^ 


+  Pj^fj^ja"'*‘^'^^'^^(loga)'^'*’^]  d|d2[ 


a^^(loga)'^  • 


Here  we  observe  that  indices  arc  related  as, 

n+l+i+j  -  p+q,  and  J  =  q» 

n+2+i+J  «•  p+q,  and  J  +  1  “  q, 


therefore  in  both  cases  n  p-l-i  so  that  we  find, 


(IV.6) 


"Pq  ^  -p-l-l-l  q  -p-l-l‘i,q-l 

V/riting  the  right  hand  side  of  (III. 6)  in  a  power  series. 


h(5!,a»a)  -  h  (x,z)  a  , 
n-0  " 

we  can  express  the  integral  equation  in  power  series  as, 

(IV.7)  “"S- 

^  ^ 5=0  po^ 


If  q  -  0,  (IV.7)  yields. 


<"'•«)  "po  -  "p  -  S  Vl-Ao  > 


and  if  q  /  0, 


(IV. 9) 


^pq  “  ’  2^  //s  ^  ^®p-l-i^iq'*’^p-l-i*’iq-l^  * 


These  expressions  are  recursion  formulas  which  permit  the  determination  of 
the  coefficients  fj^j(3t,z)  in  (IV.4)  by  means  of  iteration.  For  reference, 
we  write  down  the  first  few  terms,  from  (IV, 8) 


^oo  "  ^o  * 


''lo  '  H  ~  “o^^oo  ‘’^“5  ■  ’’l  ■  2ii//s  ?  ’ 

(IV.IO)  -  hj  -  didi 

-  "2  -  islls  (“"-iT-ldfi  f)"  ♦  ?  ‘"l-  jj/Zs  ?  «  ”  «?<'?  • 

<■30  ■  S  -  -slli  ‘"2^00*  “o''2o)  "I'’? 
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-  h  -  5?  fx-x-  S  //?  7  ) 

♦  ?{ V  I^V  ^"x 

We  find  from  (IV. 9)  f^q  ■  0  for  all  q,  and 

hi---s  Hi  ‘“.'ox*  V«.>  <«''?  ■  ^// r  ' 

(IV.XX)  fij  -  -  ^  //;  P/„i)  «  -  0. 

fjX  -  -  S  //S  ‘“X^OX*  VXX*  »ofxo>  <*r<*2 

-  -^//st  7<  57//s  V?-*5  )-<v  ^//,  ))  d?<.i . 

The  above  process  ultimately  leads  to  (IV. 4).  V;e  say  f(5t,zta)  has  an  es¬ 
timate  of  degree  (i«J)  lf» 

f(X,Sta)  -  P(a,  a  loga)  +  o[a^(a  loga)*^], 

where  P  is  a  polynomlnal  of  degree  (i,j).  From  (IV. 3)  we  have, 

aG  -  -2  a^loga  +  o(a^loga), 

hence  we  see  that  the  product  of  aC  with  a  polynomlnal  of  degree  (l,J) 

is  a  polynomlnal  of  degree  (i+l,J+l)  plus  terms  o[a^^^(a  loga)^*^]. 

Suppose  then  that  we  have  shown  f(x,z>a)  to  have  an  estimate  of  degree 

(i,j).  Substituting  this  estimate  in  the  integral  in  (IV.7)  we  obtain 

a  polynomlnal  of  degree  (i+l,J+l)  plus  terms  of  o[a^^^(a  loga)^^^],  with 

coefficients  determined  by  the  known  quantities  f  ...,  f. ..  Right 

OO,  IJ 

hand  side  of  (IV.7)  has  estimates  of  degree  (i,0)  for  all  i,  hence  (IV. 7) 
yields  for  f(x,z:a)  an  estimate  of  degree  (i+l,J+l).  Then  substituting 
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the  series  (IV./+)  and  (IV, 3)  in  the  inte;^*al  representation  (Ill.l)  it 
can  be  seen  that  V(5r,y,S:a)  has  the  form, 

aV(J[,?,55:a)<v9  !>  V..(Sr.y.z)  a^(a  loga)'^  . 

i-0  j-0  ^ 

Now  retaining  the  terms  with  h^  in  the  coefficients  f^j(x,z) 
we  write, 

(IV. 12)  f(3f,z:a)co  f^^+  loga)f22^+ 

'*«*•«  ^oo  ”  ^o  * 

^lo  «  ■  2^//?  /  » 

(IV.13)  f2^  *  -  C(in.^.log  f)h^  -  5^/4  J°d|d?[  )]d?dE  , 

%  •  ^//s  ^2^/4  ^  *  2^/4  f  » 

^30  *  *  2^/4  f)  ^74  F°  d|dt 

*  ^^/4  §>V  ^/4r  • 

According  to  (11,12)  h^  is  equal  to  1  for  the  heave,  z  for  the  roll,  and 
X  for  the  pitch.  Hence  the  integral  Jf h^d|[d^  becomes  the  water  plane 
area  in  the  case  of  heave,  and  vanishes  in  other  cases. 

Substituting  the  asymptotic  solution  (IV, 12)  in  (111,3/4)  "  (111.36) 
we  obtain  the  first  non-vanishing  term  of  the  normalized  added  mass, 
added  moment  of  inertia,  and  damping  factor  as, 

^  -  ‘  a/4  ^  ^/4  y  'I?  )  d^cdz  -  0(1), 

%  “  "  a/4  ^  §“/4  ^  -  a  S^  =  0(a), 
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(iv.u)  3 

"x  -  - 1//»  t  lr//j  <  I//5 1  '‘5''^  -  «(*“>• 

I*  -  ■  i//j  <  2j//i  ¥  "y?  >  * 

"x  ”  r//s  t  ill's  ‘  ills  7  1 « <^<15  ■ 

Note  that  N  is  equal,  to  the  product  of  the  parameter  a  and  the  square 

y 

of  normalized  water  plane  area  S.  As  the  parameter  a  tends  to  zero,  N  , 

w 

H  ,  and  H  vanish  while  K  ,  I  ,  and  I  approach  non-zero  constants. 

X  2  y  X  2 

This  is  to  be  compared  wj  th  the  situation  in  two  dimensional  case  where  N 

V 

becomes  a  constant  while  K  tends  to  infinity  as  shown  in  [5]  and  [6]. 

y 

These  estimates  Indicate  that  the  strip  method  shown  in  [5]  is  not  expect¬ 
ed  to  yield  accurate  results  for  low  frequency. 


V.  Numerical  Procedure 

V/e  suppose  the  surface  of  a  flat  ship  to  be  an  ellipse  piven  by 
(x/a)^+  (z/b)^  »  1,  or  in  normalized  co-ordinates  3^  “  1  with 

b  ■  b/a  ,  and  develop  a  numerical  procedure  by  which  the  value  of  un¬ 
known  function  fCIT,”?)  in  the  intep:ral  equation  (III .6)  can  be  determined 
approximately.  V/c  replace  the  ecjuation  (III. 6)  by  a  set  of  linear  equa¬ 
tions  relatinp  the  values  of  f(x,'z)  at  chosen  pivotal  points  on  the 
elliptic  surface.  Then  the  surface  integrals  in  the  linear  equations  are 
evaluated  by  r/impson's  rule  /'iven  by, 

/d  (h  /d  n  q  P 

(V.l)  I  f(x,z)dxdz«t  [h  C,f(x  ,z)]dzai  k  ^  C^f ( 1 

'  c  '  /  c  1  ■‘0  m-'^  1“0 

for  a  continuous  function  f(x,z).  Here  h  “(b-a)/p,  k-(d-c)/q,  and 

C3  takes  the  values  I/3  for  1=0  or  p,  U/'i  for  l--2n,  and  2/3  for  l=2n+l, 
respectively.  D  takes  the  values  I/3  for  m  O  or  q,  /*/ 1  for  m-2n,  .and 
r/3  for  m  2nfl,  respectively. 

.Substitution  of  (111.16)  ir  (111.6)  yields. 


28 


2[  f  j.(  3r,?) +if j  ( ST,?)  ]+  //g[  f ^(|,?)+if ^ (IX )  ]{f  -iia[ Y^(  aP)+S^(  a?) -i2J^( a?)  ]} 

-  2h(3?,S)  on  y  -  S(x,5). 

If  we  write  the  real  and  lmaginai*y  parts  separately,  the  following  pair  of 
equations  will  result. 

2f^(x,z)+  f  “a  log(a?)“  |^(aF)]-fj^(|,?)naJ^(ar)jd|d^  -2h(x,z), 

(V.2) 

2fj^(3f,z)+  f/j(j{fp(?»C)"aJo(ar)+f^(|,^)[  ^  -a  log(aF)-  |^(a?)]|d^d^  -  0, 
where  R(ar)  -■  Y  (a?)  +  S  (ar)  -  ~  log(ar). 

O  On 

We  now  establish  a  lattice  on  the  elliptic  surface  S  by  dividing 
the  long  axis  into  eight  equal  intervals  h,  and  the  vertical  ordinates 
parallel  to  the  short  axis  into  four  equal  intervals  k(3r),  that  is, 
h  «■  lA  and  k(x)  -  b  /l-x'  /2.  In  the  course  of  additional  computations 
for  improving  numerical  results  these  intervals  are  bisected  to  yield  a 
finer  grid.  However  we  only  present  here  the  procedure  for  original  la¬ 
ttice,  In  this  instance,  each  pivotal  point  can  be  identified  by  the  co¬ 
ordinates  Xj  and  Zj  where  Xj^-(i-4)A  i“0,l,...8,  and 

J-0,l,...iV.  To  determine  the  values  of  f(x,z),  we  consider 
the  equation  (V.2)  only  at  thirteen  pivotal  points  contained  in  one  qua- 
dreint.  [Note  that  for  the  fine  lattice,  the  pivotal  points  contained  in 
one  quadrant  are  rrrty-one.]  V/hen  the  values  of  f(Sf,z)  at  these  points 
are  known,  the  symmetry  or  anti-symmetry  properties  of  the  function  will 
enable  us  to  determine  the  values  at  the  rest  of  pivotal  points. 

As  F  denotes  the  distance  from  a  fixed  pivotal  point  (x^,Zj)  to 
any  pivotal  point  (5»?)  ■  where  1=0,1, ...8,  and  m=0,l,..,^, 

we  find  the  integrands. 
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and  l0RaJ(x^-|)2+(2j-^" 

possess  a  sinpularity  at  (Xj^,Zj)-(Xj^,z^).  Nevertheless,  the  singular 
Integrals  associated  vdth  these  integrands  do  exist.  We  vdll  presently 
show  how  these  integrals  can  be  evaluated. 

(i)  Treatment  of  the  integral  ^ 

Il(Xi*«j)-// - d^d^  . 

We  choose  the  shortest  distance  r^  between  any  two  neighboring 
pivotal  points  and  draw  a  circular  region  about  the  fixed  pivotal  point. 
If  the  integral  Ij^(x^,Zj)  i-  evaluated  over  the  region  S-6,  and  the  ex¬ 
cluded  circular  region  6  separately. 


(V.3) 


whore  the  function  f(3r,?)  is  regarded  as  a  constant  over  the  region  6, 
Since  the  integral  over  the  circular  region  becomes. 


nr  , 

°  i  d?  rd(3 
o 

wc  ODiain, 

(V,4)  Ij^(5!'j,Zj)  x  2nr^f(xj,Zj) 


Now  suppose  the  integral  Ij^(Xj^,Zj)  is  evaluated  by  (V.l)  assigning  a 
fictitious  value  r^  for  r  at  the  singularity. 
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where  the  notation  >q  denotes  the  double  eunnatlon  with  9 

inside  the  circle.  C.  and  D  are  the  coefficients  of  Simpson*  s  rule  in 

1  m 

the  X-  and  z-directions,  respectively. 

Hence  we  find> 


(V.5)  ll  ■  .  dld^g  nrpf(SEj^,gj) 


O 


'=l''D.k(*i> 


«*1‘V 


Substitution  of  (V.5)  in  (V.4)  yields, 

(V.6) 


8  L 


For  pivotal  points  on  the  boundary,  the  correction  term  «r^f(^,Zj)  will 

be  approximated  by  one  half  because  the  circular  region  drawn  with  the 

radius  r  about  such  a  pivotal  point  does  not  make  a  full  circle  as  the 
o 

region  will  be  sliced  off  by  the  boundary  of  the  ellipse. 

(11)  Treatment  of  the  Integral 

logay(xj-|)*+(zj-'5)*  . 

If  the  integral  evaluated  separately  over  the  region 

S>6  and  6, 

(V.7)  ^  y(x^-|)M*j-5)*  d^dt 

♦  f(xj^,zJ  /4  loga  j(x^-|)*+(»4-?)®  d^dJJ  , 


where  f(S,z)  is  again  regarded  as  constant  over  the  region  6. 
Sir.ee  the  integral  over  the  circular  region  becomes, 
f2n.tr^ 

I  log  a?  dr  rde  -  nr^dogar^  ”  2 

*  o  '  o 
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we  have, 


(V.8)  l2(Xj^,Zj)»nrJ(logar^- f (l,?)loga  j(x^-?)®+( Zj-?)^d^dt. 

Assigning  a  fictitious  value  for  ?  at  the  singularity,  and  evaluating 
l2(j^«2j)  by  (V.l)  we  obtain, 

//s-A  d|d2[  +f(x^,Zj)logar^  //^  dfd^ 

Therefore,  ____________ 

(V.9)  f(|,^)loga  d|d^  »  -nrJlogar^f(x^,Zj) 


1^ 


C j^hD^k  ( 5^  )  loga  ^  (  Sr^  )  “+ ( z  -a^ )  ^  fCxj^,?^^^). 


Substitution  of  (V.9)  in  (V.8)  yields, 


(V.IO)  l2(Xi,?j)  *  -  §**0 


C^hD^k(Srj)loga  A^i'V^(*J*^m)" 


1*0 


Now  the  application  of  Simpson's  rule  (V.l)  together  with  the 
singular  integral  formulae  (V.6)  and  (V.IO)  enable  us  to  reduce  the  set 
of  integral  equation  (V.2)  to  twenty-six  (for  the  original  lattice  or 
eighty-two  for  the  fine  lattice)  linear  equations  relating  the  values  of 
f^(Sr,z)  and  f^(x,z)  at  chosen  pivotal  points. 

For  convenience  we  use  the  following  notations, 

f^J  =  f(xj,Zj),  f^-  f(x^,zj, 

(V.ll) 

K,  C,h  D  k(x,),  and  h^"^ 
im  1  m  1 


-  MSTj.Tj). 


We  write  (V,2)  as. 
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t2  *nar„(l-  |£o)lf«  *  f  Kj„[ 


1-0  m-0 


(V.12) 


+nar^(l- 


where  n  is  a  function  of  (i,j)  which  takes  the  value  of  1/2  or  1  .  Note 

that  amon^  the  coefficients  of  Simpson's  rule  the  relation,  K,  “K„  , 

X  fiti  7“X  ,in 

■K,  .  _"K-  I  e__  1  holds.  We  further  write, 
li5-m  9-1, 5"m  ’  ’ 


(V.13) 


C»ar^(l  - 


Next  we  shall  investif»atc  the  charecteristics  of  linear  equations  associat¬ 
ed  with  the  heave,  roll  and  pitch,  respectively! 

(1)  The  case  of  heave. 

Since  the  function  f(5l[,z)  satisfies  the  symmetry  relations, 

_l,m  ,9-l,in  .  ,l,m  _l,5*fli  -  u 

f  *  -  f  *  ,  and  f  *  -  f  *  for  heave. 


the  double  sums  appearing  in  (V.12)  can  be  expressed  as, 

(V.U)  X  ^  K,  H^f^ 

1-0  m-0 


1-0  m-0 


From  (V.13)  and  (V.14)  the  linear  equations  (V.12)  for  heave  becomes. 
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Jl  2 


(2  +nC)fj;'^  + 

1^0  in“0 


Iffl 


(V.15) 


1-0  m-0 


(2  +nC)fJ'^  + 


Im 


•a 


1-0  m-O 


The  coefficient  matrix  of  (V.15)  Is  antl-symnetrlc ,  and  Its  elements  In 
the  first  and  the  fourteenth  colvunns  and  row  will  vanish  except  the  main 
diagonal  elements  because  In  the  process  of  the  double  Integration  by 
(V.l)  the  function  at  the  both  ends  of  the  major  axis  are  not  taken  Into 
account. 

(11)  The  case  of  roll. 

For  roll  the  function  f(52',?)  satisfies  the  symmetry  relation, 

-l.m  ,9*1, m  ...  .  .  .  ...  -l.m  -9*l,m  ,l,5*in  -9»1»5*« 

f  *  -  f  *  ,  and  the  anti-symmetry  relation  f  '  -f  *  "-f  '  — f  *  , 

therefore  the  double  suns  appearing  In  (V.12)  for  the  roll  becomes, 

(V.16) 


1-0  m-0  1-0  m-0 

From  (V.13)  and  (V.l6)  the  linear  equations  (V.12)  for  the  roll  becomes, 

(2  +nC)f^J 

1"  ^  ^  A  f\  Am  A  A  A  A  r 

1-0  m-0 


1-0  m-0 


2  +nC)fjJ  + 


1-0  m-0 


(V.17) 


1*0  m-0 

Here  inspection  of  (V.l6}  shows  that, 


L  2  , 

-a"  ^21  -  0. 

tTSs:  zny:  m  ^  2  2  2  r 


(V.18) 


^  for  m  -  2 


1«0  wO 


Therefore  the  coefficient  matrix  of  (V,17)  is  also  anti-'Synmeti*ic  and 
in  addition  to  vanishing  first  and  fourteenth  columns  and  rows,  from 
(V«18)  it  has  zeroes  in  the  (4'*'3n)th  and  (17‘*-3in)t>h  columns  and  rows 
except  the  main  diagonal  elements,  m  here  assumes  the  values  0,1,2,3< 
(i'.i)  The  case  of  pitch. 

For  pitch  the  function  f(3t,z)  has  the  symmetry  property,  f^*"™ 


■f^*^^,  and  the  anti-symmetry  property,  f^*'"-f^*^^— f*^"^*"*»-f^  1»5*«^ 
therefore  the  double  sums  appearing  in  (V,12)  for  the  pitch  becomes. 


(V.19)  •^l,n(H^*"‘+H^*5in^9-l,m^9-l,5-m)jlm^ 

1«0  m-0  1-0  m-O 

From  (V.I3)  and  (V.19)  we  write  the  linear  equations  (V.12)  for  pitch  as. 


(2  +nC)f^J  +  J^SKij„(Hi*"'+Hj»5in.H9-l,m^9-l,5-in)j.L 

1-0  m-0 


Im 


-a*  Kj^(Hj»'"+Hj»5-m^9-l,m^9-l,5-mjj^..2x, 


1-0  m-O 


(V.20) 


(2  +nC)fJ'^  +  To 


IrO  m-O 


1-0  m-0 


Here  it  can  be  seen  that. 


(v.?l) 


_  q  for  1  -  4, 


(Hl»m+Hl*5-fn_H9-l,m^9-l,5-m)_  ^  for  i  -  4. 
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m-O 


Therefore  the  coefficient  matrix  of  (V.20)  is  anti-8ynmetric>  and  from 
(V.2I)  it  hae  vanishing  elements  in  the  (ll-Hn}th  and  (24'Hn)  th  columns 
and  rows  except  the  main  diagonal  elements,  m  here  assumes  the  values 
0,1, 2« 

The  integral  equation  (V.6)  describing  the  forced  oscillation  of 
a  flat  ship  is  thus  replaced  by  three  sets  of  linear  equations*  For  vari* 
ous  values  of  the  parameter  a*  (^^a/g  ^2n&/K  ,  we  can  determine  the  values 
of  f(St,t)  by  solving  these  equations*  Actual  steps  of  the  computation  work 
for  obtaining  the  solution  consist  of  t 

(1)  Determination  of  the  distance  r  from  individual  pivotal  point  to  any 
pivotal  point  on  the  given  lattice* 

(2)  Calculation  of  functions  appearing  in  the  coefficients  of  the  linear 
equations  either  by  direct  evaluation  or  by  interpolation  from  the  given 
table* 

(3)  Summation  of  the  coefficients  and  grouping  of  the  matrix  in  accordance 
with  the  type  of  motion* 

(4)  Numerical  solution  of  linear  equatlceis  either  by  the  elimination  process 
or  iteration  process* 

In  step  (1)  when  the  distance  F  becomes  zero,  it  is  replaced  by  the 
smallest  distance  between  two  neighboring  pivotal  points  r^  ,  and  in 
step  (2)  the  functions,  l/r  and  logaF  are  evaluated  directly  and  the  fUnc- 
tlon  R(a?)-  y  (a?)  +  S^(ar)  -  —  logaF  and  J  (a?)  are  evaluated  by  means  of 
parabolic  interpolation  from  a  pre-arranged  table  [Table-1]*  In  step  (4) 
successive  elimination  of  unknowns  based  on  the  algorithm  of  Gauss  is  used. 
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Table  1 


Input  for  Determination  of  Coefficients  of  Linear  Equation 


a? 

Ri&lf) 

Jo(ar) 

af 

R(ar) 

Jo(a?) 

ar 

R(8?) 

Jo(a7) 

0,00 

-0.0738 

1.0000 

2.25 

0.773U 

0.0828 

I4.50 

-1.2107 

-0.3205 

0.05 

-O.OliOl 

0.999U 

2.30 

0.7552 

0.0555 

U.55 

-1.2U87 

-0.3087 

O.IO 

-O,oob5 

0.9975 

2.35 

0.7303 

0.0288 

I4.6O 

-1,2851 

-0.2961 

0.15 

0.0320 

O.99I4I4 

2.U0 

0.7036 

0.0025 

U.65 

-1.3200 

-0,2830 

0.20 

0.070b 

0,9900 

2.U5 

0.6751 

-0.0232 

U.70 

-1.3533 

-0.2693 

0.25 

0.1085 

O.98I4I4 

2.50 

0.6l4it8 

-0.0(4814 

U.75 

-1.3850 

-0.2551 

0.30 

0.1ii86 

0.9776 

2.55 

O.612I4 

-0.0729 

li.30 

-I.U151 

-O.2UOU 

0.35 

0.1681 

0.9696 

2.60 

0.578U 

-0.0968 

U.35 

-1.UU35 

-0,2253 

0.1i0 

0.227li 

O.960I4 

2.65 

0.5I428 

-0.1200 

U.90 

-1.)4702 

-0.2097 

o.)45 

0.267U 

0.9500 

2.70 

O.505I4 

-O.II42I4 

U.95 

-1.U951 

-0.1938 

0.50 

0.3059 

0.9385 

2.75 

0.14667 

-O.I6I4I 

5.00 

-1.5183 

-0.1776 

0.55 

0.3551 

0.9258 

2.80 

0.)i263 

-0.1850 

5.05 

-1.5398 

-0,1611 

0.60 

0.3831 

0.9120 

2.85 

O.38I46 

-0.2051 

5.10 

-1.559U 

-0.IUU3 

0.65 

O.I4209 

0.8971 

2.90 

0.3U15 

-0.22h3 

5.15 

-1.5772 

-0.127U 

0.70 

O.U579 

0.3812 

2.95 

0.2972 

-0,21426 

5.20 

-1.5933 

-0.1103 

0.75 

O.I49UU 

O.86I42 

3.00 

0.2518 

-0.2601 

5.25 

-1.6075 

-0.0931 

0.30 

0.5293 

0.3U63 

3.05 

0.2052 

-0.2765 

5.30 

-1.6199 

-0.0758 

0.85 

0.5631 

O.327I4 

3.10 

0.1576 

-0,2921 

5.35 

-1.6306 

-0.0585 

0,90 

0.5957 

0.8075 

3.15 

0.1093 

-0.3066 

5.U0 

-1.639U 

-0.0(412 

0.95 

0,6271 

0.7868 

3.20 

0.0600 

-0.3202 

5.U5 

-I.6U63 

-0.02U0 

1.00 

0.6570 

0.7652 

3.25 

0.0099 

-0,3328 

5.50 

-1.6515 

-0,0068 

1.05 

0,6852 

0.7U28 

3.30 

-0,0507 

-0.3UU3 

5.55 

-1.6552 

0.0102 

1.10 

0.7121 

0.7196 

3.35 

-0.0917 

-0.351*8 

5.60 

-1.6569 

0.0270 

1.15 

0.7372 

0,6957 

3.U0 

-O.II435 

-O.36I43 

5.65 

-1.6570 

0.0U36 

1.20 

0.7606 

0.6711 

3.U5 

-0.1955 

-0.3727 

5.70 

-1.6555 

0.0599 

1.25 

0.7822 

O.6I459 

3.50 

-0.2I478 

-0,3801 

5.75 

-1.6523 

0.0760 

1.30 

0.8019 

0.6201 

3.55 

-0,3001 

-0.3865 

5.80 

-1.6U75 

0,0917 

1.35 

0.8197 

0.5937 

3.60 

-0.3526 

-0.3918 

5.35 

-1.6U12 

0.1071 

1.U0 

0.835U 

0.5669 

3.65 

-0,)4052 

-0.3960 

5.90 

-1.633U 

0.1220 

1.U5 

O.8I492 

0.5395 

3.70 

-0.I457U 

-0.3992 

5.95 

-1.62UO 

0.1366 

1.50 

0.8610 

0.5118 

3.75 

-0,5096 

-0.l401i4 

6.00 

-1.6135 

0.1506 

1.55 

0.8706 

0.14838 

3.80 

-0.5616 

-O.I4O26 

6.05 

-1.601U 

0.16U2 

1.60 

0.3782 

O.U5514 

3.85 

-0,6131 

-0.14027 

6.10 

-1.5880 

0.1773 

1.65 

0,8836 

O.I4268 

3.90 

-O.66I4I 

-0,14010 

6.15 

-1.5735 

0,1898 

1.70 

0.8868 

0.3980 

3.95 

-O.73I47 

-O.I4OOO 

6.20 

-1.56U1 

0.2017 

1.75 

0.8879 

0.3690 

I4.00 

-O.76U4 

-0.3971 

6.25 

-1.5U09 

0.2131 

1.80 

0,8867 

O.3I4OO 

I4.05 

-0.8137 

-O.393I4 

6.30 

-1.5230 

0.2230 

1.85 

0.8833 

0.3109 

U.io 

-0.8620 

-0.3887 

6.35 

-I.50UI 

0.2339 

1.90 

0.3777 

0.2818 

U.15 

-0.9093 

-0,3831 

6.U0 

-I.I48U2 

0.2U33 

1.95 

0.8700 

0.2528 

U.20 

-0.9559 

-0.3766 

6.I45 

-I.U637 

0,2521 

2,00 

0.8601 

0.2239 

I4.25 

-l.OOlli 

-0.3692 

6.50 

-1.UU21 

0,2601 

2.05 

O.8I476 

0.1951 

I4.30 

-I.OI457 

-0.3610 

6.55 

-1.U200 

0.267U 

2.10 

0.8335 

0.1666 

U.35 

-1.0091 

-0.3520 

6,60 

-1.3973 

0.27U0 

2.15 

0.8170 

0.1383 

b.iiO 

-1.1308 

-O.3I423 

6.65 

-1.3730 

0.2799 

2.20 

O.798I4 

O.IIOI4 

l4,l45 

-1.1715 

-0.3318 

6.70 

-1.3501 

0.2351 
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Table  1 

Input  for  Determination  of  Coefficients  of  Linear  Equation 


er 

R(a?) 

ar 

R(ar) 

ar 

R(ar) 

Jo(a?) 

7.55 

-0.9U37 

0.2593 

9.80 

-1.1787 

-0.2323 

12,05 

-1.9763 

0,0588 

7.60 

-0.9238 

0.2516 

9.85 

-1.2071 

-0.2366 

12.10 

-1.9719 

0.0697 

7.65 

-O.90I48 

0.2li3U 

9.90 

-1.2350 

-O.2I4O3 

12.15 

-1.9663 

0,0803 

7.70 

-0.8865 

0.23U6 

9.95 

-1.2631 

-0,2U3U 

12.20 

-I.95I48 

0.0908 

7.75 

-0.8691 

0.2252 

10.00 

-1.2915 

-0.2159 

12.25 

-1.9520 

0.1009 

7.80 

-0.8528 

0.215U 

10.05 

-1.3199 

-0.2)478 

12.30 

-1.9U35 

0.1108 

7.85 

-0.8375 

0.2051 

10.10 

-1.3li8U 

-0.2I490 

12.35 

-1.9339 

0.1203 

7.90 

-0.8233 

O.I9U4 

10.15 

-1.3770 

-O.2I496 

I2.I4O 

-1.923)4 

0,1296 

7.95 

-0.8103 

0.1832 

10.20 

-I.U055 

-0,2li96 

I2.I45 

-1.9121 

O.138I4 

8.00 

-0.7983 

0.1717 

10.25 

-1.)4338 

-0,2)490 

12.50 

-1.8997 

O.II469 

8.05 

-0.7876 

0.1597 

10.30 

-I.I4620 

-0.2)477 

12.55 

-1.8666 

0.1550 

8.10 

-0.7780 

0.11i75 

10.35 

-1.14899 

-0,2)458 

12.60 

-1.8730 

0.1626 

8.15 

-0.7698 

0.1350 

10.U0 

-1.5176 

-O.2I43I4 

12,65 

-1.858)4 

0.1698 

8.20 

-0.7627 

0.1222 

10.U5 

-1.5U49 

-O.2I4O3 

12.70 

-1.8)433 

0.1766 

8.25 

-0.7570 

0.1092 

10.50 

-1.5718 

-0,2366 

12.75 

-I.827I4 

0,1829 

8.30 

-0.7526 

0.0960 

10.55 

-1.5985 

-O.232I4 

12.80 

-1.8110 

0.1887 

8.35 

-0.7U9U 

0.0826 

10.60 

-1.62)42 

-0.2276 

12.85 

-1.79)40 

0.19)40 

8.U0 

-0.7U77 

0.0692 

10.65 

-1.6U96 

-0.2223 

12,90 

-1.7767 

0.1988 

8Ji5 

-0.7U72 

0.0556 

10.70 

-1.67U3 

-0.216)4 

12.95 

-1.7588 

0.2031 

8.50 

-0.7U80 

O.OliW 

10.75 

-1.6983 

-0.2101 

13.00 

-I.7I4O6 

0.2069 

8.55 

-0.7U9U 

0.0283 

10.30 

-1.7216 

-0.2032 

13.05 

-1,7221 

0.2102 

8.60 

-0.7538 

O.OII16 

10.85 

-I.7I4I42 

-0.1959 

13.10 

-I.703I4 

0.2129 

8.65 

-0.7586 

0.0010 

10.90 

-1.7660 

-0.1881 

13.15 

-1,68)43 

0.2151 

8.70 

-0.761»9 

-0.0125 

10.95 

-1.7868 

-0.1798 

13.20 

-1.6652 

0.2167 

8.75 

-0.7725 

-0.0259 

11.00 

-1.8068 

-0.1712 

13.25 

-1.6U58 

0.2178 

8.80 

-0.7812 

-0.0392 

11.05 

-1.8259 

-0.1622 

13.30 

-1.6267 

0.2183 

8.35 

-O.79IU 

-0.0523 

u.io 

-I.8UI0 

-0,1528 

13.35 

-1.6118 

0.2183 

8.90 

-0.8026 

-0.0653 

11.15 

-1.8611 

-0,11430 

13.I4O 

-1,5880 

0.2177 

8.95 

-0.8l5i 

-0.0779 

11.20 

-1.8770 

-0,1330 

13. U5 

-1.5689 

0.2166 

9.00 

-0.3290 

-0.0903 

U.25 

-1.8922 

-0.1227 

13.50 

-I.5I498 

0.2150 

9.05 

-O.8I439 

-0.102U 

U.30 

-1.9060 

-0.1121 

13.55 

-1.5311 

0.2128 

9.10 

-0.8600 

-O.IU42 

U.35 

-1.9188 

-0.1012 

13.60 

-1.5126 

0.2101 

9.15 

-0.8771 

-0.1257 

ll.liO 

-1.9305 

-0.0902 

13.65 

-l,i49U5 

0,2069 

9.20 

-0.895U 

-0.1367 

11.U5 

-1.9U09 

-0.0790 

13.70 

-1.)4767 

0.2032 

9.25 

-0.911i5 

-0.11t7U 

11.50 

-1.9503 

-0.0677 

13.75 

-1.U593 

0.1990 

9.30 

-0.93U8 

-0.1577 

11.55 

-1.9586 

-0.0562 

13.80 

-1.UU25 

0.19)i3 

9.35 

-0.9560 

-0.167U 

11.60 

-1.9657 

-O.OI4I46 

13.85 

-I.I4260 

0.1892 

9.'40 

-0.9780 

-0.1768 

11.65 

-1.9713 

-0.0330 

13.90 

-l.)4l02 

0.1836 

9.U5 

-1.0008 

-0.1856 

11.70 

-1.9761 

-0.0213 

13.95 

-1.3951 

0.1775 

9.50 

-1.02U5 

-0.1939 

11.75 

-1.9806 

-0,0097 

II4.00 

-1.3805 

0.1711 

9.55 

-l.OliSe 

-0.2017 

11.80 

-1.9820 

0,0020 

II4.05 

-1.3667 

0.16)42 

9.60 

-1.0739 

-0.2090 

11.85 

-1.9826 

0.0135 

II4.I0 

-I.353I4 

0.1570 

9.65 

-1.0995 

-0.2157 

11.90 

-1.9831 

0,0250 

1)4.15 

-I.3I4I2 

0.1)493 

9.70 

-1.1257 

-0.2218 

11.95 

-1.9821 

0,036)4 

1)4.25 

-1.3187 

0.1331 

9.75 

-1.1525 

-0.2273 

12.00 

-1.9796 

O.OI477 

1)4.30 

-1.3087 

0.12)45 
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VI .  Dlscuaslon  of  Numerical  Heaults 

The  computations  are  in  terms  of  two  parameters.  The  first  is 
the  ratio  of  the  short  axis  to  the  lonf;  axis  of  the  ellipse  b  "b/a, 
and  the  second  parameter  is  the  ratio  of  the  half  len/;th  of  the  ship 
to  the  wave  lenf^th  a  -  2k&/K"  ^*a/g.  The  values  of  b  chosen  for  the 
investigation  were  1/8,  l/U,  and  1  so  that  they  represent  a  slim  ellip¬ 
se  and  a  circle  as  limiting  cases.  The  parameter  a  assumes  the  values 
n/6,  n/5,  n/i»,  n/3,  2n/5,  ii/2,  2n/3,  and  n. 

It  was  shown  in  section  III  that  for  large  argument  of  af^,  the 
asymptotic  form  of  the  OreenV  function  contains  a  trigonometric  func¬ 
tion,  that  is,  G(x,z,^,^)s  2/r  +  i2^ 2na/i^  ^/U) ^  Hence, 

the  kernel  of  the  integral  equation  fluctuates  as  the  frequency  of  os¬ 
cillation  Increases.  This  implies  that  to  assure  equal  accuracy  we  must 
take  the  grid  spacing  inversely  proportional  to  the  frequency.  For  this 
reason,  the  computations  are  carried  out  at  first  with  the  original  la¬ 
ttice  and  then  with  the  fine  lattice  which  has  the  bisected  grid  ipa- 
cings. 

For  each  combination  of  b  and  a,  three  set  of  linear  equations 
(V,15),  (V,17),  and  (V.20)  were  solved  in  order  to  determine  the  values 
of  f^(3!,z)  and  fj^(x,z)  at  chosen  pivotal  points.  Then  using  Simpson's 


wore  evaluated  by  (III .34).  SimilaiT-,  for  roll  the  normalized  added  mo¬ 
ment  of  inertia  I^  and  damping  factor  were  evaluated  by  (111.35), 

while  for  pitch  the  normalized  added  moment  of  inertia  I  and  damping 

z 

factor  H  were  evaluated  by  (III. 36).  In  Table  2  the  values  of  M  ,  N  , 

z  ^  y  y 

I  ,  H  ,  I  ,  and  H  for  a  circular  disk,  b-1,  corresponding  to  various 

X  X  T,  7. 

values  of  the  parameter  a  are  tabulated.  In  Tibles  3  and  4  the  values 

of  M  ,  N  ,  I  ,  H  ,  I  ,  and  H  for  elliptic  disks  of  the  axes  ratio  I/4, 
y  y  X  X  z  z 

and  1/8,  depending  upon  the  parameter  a  are  presented.  In  those  tables, 
the  values  within  the  parenthesis  denote  the  results  obtained  by  the 
use  of  the  original  lattice  and  the  other  values  by  the  use  of  the  fine 
lattice.  Note  that  for  the  case  of  elliptic  disks  the  results  obtained 
by  the  use  of  the  original  lattice  are  not  much  different  from  those  by 
the  use  of  the  fine  lattice.  However  in  Table  2  it  can  be  seen  that  the 


Table  2 
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Added  Mase,  Added  Moment  of  Inertia,  and  Damping  Factors 

for  Circular  Disk  ^1. 


N, 

I 

H 

I 

H 

a  -  «/6 

y 

y 

X 

X 

s 

z 

2.2U2 

0.990 

0.351 

0.032 

0.335 

0.C30 

(2.m) 

(0.977) 

(0.262) 

(0.029) 

(0.194) 

(o.o  3) 

a  -  «/5 

2.129 

1.016 

0.415 

O.O64 

0.395 

0.060 

(2.082) 

(1.006) 

(0.306) 

(0.055) 

(0.232) 

(0.044) 

a  -  n/k 

1.987 

1.026 

0.536 

0.161 

0.509 

0.150 

(1.943) 

(1.021) 

(0.385) 

(0.130) 

(0.299) 

(0.101) 

9> 

1 

;« 

> 

1.809 

1.002 

0.677 

0.651 

0.650 

0.599 

(1.774) 

(1.007) 

(0.483) 

(0.458) 

(0.387) 

(0.338) 

a  -2n/5 

1.706 

0.968 

0.126 

1.136 

0.163 

1.078 

(1.659) 

(0.970) 

(0.069) 

(0.874) 

(0.125) 

(0.683) 

a  “  k/2 

1.593 

0.910 

-0.497 

0.602 

-0.469 

0.601 

(1.582) 

(0.923) 

(-0.366) 

(0.575) 

(-0.277) 

(0.543) 

a  -2r/3 

1.473 

0.811 

-0.381 

0.167 

-0.370 

0.168 

(1.489) 

(0.809) 

(•0.336) 

(0.147) 

(-0.316) 

(0.178) 

a  -  « 

1.3a 

0.628 

-0.127 

0.031 

-0.138 

0.009 

(1.342) 

(0.533) 

(-0.151) 

Table  3 

(0.001) 

(-0.190) 

(0.017) 

Added  Mass,  Added 

Moment  of  Inertia,  and  Damping 
for  Klliptic  Disk  b-l/4. 

Factors 

M 

N 

I 

H 

I 

H 

a  ■  n/6 

y 

y 

X 

X 

z 

z 

0.298 

0.099 

0.0012 

0.00001 

0.052 

0.002 

(0.294) 

(0.097) 

(0.0018) 

(0.00001) 

(0.040) 

(0.001) 

a  -  n/5 

0.288 

0.108 

0.0013 

0.00001 

0.057 

0.003 

(0.284) 

(0.106) 

(0.0018) 

(0.00001) 

(0.044) 

(0.003) 

a  -  n/U 

0.274 

0.119 

0.0014 

0.00002 

0.065 

0.007 

(0.270) 

(0.116) 

(0.0018) 

(0.00002) 

(0.051) 

(0.006) 

a  n/3 

0.251 

0.129 

0.0015 

0.00005 

0.084 

0.020 

(0.248) 

(0.127) 

(0.0018) 

(0.00005) 

(0.065) 

(0.016) 

Table  3 
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Added  Maes,  Added  Moment  of  Inertia,  and  Dampinf^  Factors 
for  Elliptic  Disk  h^l/U. 


M 

N 

I 

I 

H 

y 

y 

X  X 

t 

e 

a  -2»i/5 

0.235 

0.133 

0.0017  0.00010 

0.100 

0.043 

(0.233) 

(0.131) 

(0.0018)  (0.00008) 

(0.078) 

(0.033) 

a  -  n/2 

0.215 

0.135 

0.0019  0.00020 

0.102 

0.112 

(0.213) 

(0.133) 

(0.0019)  (0.00017) 

(0.082) 

(0.083) 

a  -2n/3 

0.187 

0.132 

0.0025  0.00058 

-0.041 

0.139 

(0.187) 

(0.131) 

(0.0021)  (0.00044) 

(-0.020) 

(0.121) 

a  “  11 

0.15A 

0.117 

0.0032  0.00312 

-0.037 

0.014 

(0.155) 

(0.118) 

(0.0023)  (0.00184) 

(-0.040) 

(0.018) 

Table  4 

Added  Mass,  Added  Moment  of  Inertia,  and  Dampin/^  Factors 

for  Elliptic  Disk  b-l/8. 

M 

N 

I  H 

I 

H 

y 

y 

X  X 

z 

z 

a  •  n/6 

0.105 

0.029 

0.00011  0.0000001 

0.020 

0.0004 

(0.107) 

(0.028) 

(0.00027)  (0.0000001) 

(0.018) 

(0.0003) 

a  ■  a/5 

0.102 

0.032 

0.00011  0.0000002 

0.021 

0.001 

(O.IOA) 

(0.031) 

(0.00027)  (0.0000002) 

(0.019) 

(0.001) 

a  -  n/U 

0.098 

0.036 

0.00011  0.0000003 

0.023 

0,002 

(0.100) 

(0.035) 

(0.00026)  (0.0000003) 

(0.022) 

(0.001) 

a  -  n/3 

0.090 

o.oa 

0.00012  0.0000008 

0.027 

0.004 

(0.092) 

(0.039) 

(0.00025)  (0.00(XX)07) 

(0.026) 

(0.003) 

a  -2a/5 

0.085 

0.0A3 

0.00013  0.0000013 

0.031 

0.008 

(0.087) 

(o.oa) 

(0.00025)  (0.0000013) 

(0.030) 

(0.007) 

a  -  n/2 

0.077 

0.045 

0.00013  0.0000027 

0.037 

0.018 

(0.079) 

(0.043) 

(0.00024)  (0.0000025) 

(0.035) 

(0.016) 

a  *211/3 

0.066 

0.041. 

0.00015  0.0000066 

0.027 

0.046 

(0.068) 

(0.042) 

(O.oa)23)  (0.00(XX)59) 

(0.027) 

(0.043) 

a  -  n 

0.051 

0.039 

0.00020  0.0000252 

-0.017 

0.023 

(0.053) 

(0.037) 

(0,00021)  (0.0000208)( -0.016) 

(0.022) 

m 

fin«  lattice  yield  much  improved  results  since  I  and  I  ,  and  H  and  H 

X  Z  X  z 

must  be  equal  in  the  case  of  a  circular  disk. 

In  Figure  1  and  Figure  2  the  dependence  of  and  on 

Ay  "  y 

the  parameter  a  are  presented.  In  Figure  3  and  Figure  U  the  quantities 
na^  Ra* 

-^I  and  H  are  plotted^as  functions  of  the  parameter  a.  Uimilary, 

Ax  AX  11^52  Ra^ 

the  quantities  — r— I_  and  ^  ti  are  plotted  in  Figure  5  and  Figure  6. 

The  multiplication  factor  ra^/A  for  the  ordinates  represents  the  ratio 
of  the  area  of  a  circle  having  the  half  length  of  the  ship  a  as  the  ra> 
dlus  to  the  area  of  the  water  plane  of  the  ship  under  consideration ,  and 
was  introduced  in  order  to  make  the  curves  comparable.  Note  that  these 
curves  are  obtained  from  the  results  of  the  fine  lattice. 

The  curves  for  the  circular  disk  b-1  in  Figure  1  and  Figure  2  compare 
very  closely  to  the  corresponding  curves  in  Figure  6  and  Figure  7  in  [5] 
which  were  obtained  by  treating  the  circular  disk  as  an  axial  symmetric 
two-dimensional  configuration. 

To  ascertain  the  accuracy  of  the  results  represented  by  the  curves  for 
S«l//»  and  b«l/8  in  Figure  2,  the  normalized  damping  factor  N'  were  com- 

y 

puted  by  the  following  formula  based  on  the  strip  method, 


(VI.l)  NJ^  -  2  (1-5*)  NgCl  a  !/l5*)d3f  , 

where  the  value  N^C^  a  /l-x*)  are  taken  from  Figure  i*.  in  [5]  using  the 
relation,  a  ,ll-X^)  ■  .  As  shewn  in  Table  5»  for  low  fre¬ 
quency  the  strip  method  does  not  yield  a  satisfactory  results.  For  this 
range  the  present  method  should  ,"ive  accurate  results. 


Table  5 


Comparison  of  Damping  Factor  of  Klliptic  Disks  for  Heave  Evaluated  by 


Integral  Iv^uation  Method  and  Strip  Method. 


For  b’»l/4 

For  b-1/8 

• 

N 

n" 

y 

y 

y 

y 

y 

«/3 

0.127 

0.129 

0.183 

(Not  Available) 

n/2 

0.1)3 

0.135 

0.166 

0.043 

0.044 

?n/3 

0.131 

0.132 

0.153 

0.0/*2 

0.044 

0.049 

0.046 


U2 

Here  in  Table  5,  N°  denotes  the  value  obtained  by  the  use  of  the  original 
f  ^ 

lattice  and  denotes  that  by  the  use  of  the  fine  lattice. 

An  additional  check  for  the  results  of  the  heave  of  a  circular 
disk  can  be  made  by  comparing  the  values  of  f^(X,S)  and  fj^(3f,2)  at  pivo¬ 
tal  points  of  the  equal  radial  distances,  V.'e  compare  these  values  at  the 
tip  and  the  half  radial  distance  on  the  3f-axio  with  those  at  the  tip  and 
at  the  half  radial  distance  on  the  ?-axis  in  Table  6  using  the  results  of 
the  fine  lattice.  At  higher  frequency  the  agreement  was  found  to  be  unsa¬ 
tisfactory  presumably  due  to  the  use  of  non-square  grid  which  is  primarily 
designed  for  the  elliptic  disk. 


Table  6 

Comparison  of  Real  and  Imaginary  Parts  of  Density  for  Heave  at  Pivotal  Points 
of  Equal  Radial  Distances  on  Circular  Disk, 


fjx, 

M) 

f^(3l 

,2) 

At 

Tip  on  5, 

Tip  on  2,1/2  on  5c, 

,1/2  on  2, 

Tip  on  St, 

,Tip  on  z. 

1/2  on  X, 

,1/2  on  z 

a 

■  n/6 

0.727 

0.712 

0.587 

0.596 

-0.185 

-0.181 

-0.160 

-0.163 

a 

■  n/5 

0,6% 

0.678 

0.526 

0.536 

-0.232 

-0.226 

-0.196 

-0.200 

a 

-  n/U 

0,655 

0.635 

0.442 

0.452 

-0.301 

-0.292 

-0.245 

-0.251 

a 

-  n/3 

0.604 

0.579 

0.312 

0.320 

-0.409 

-0.393 

-0.315 

-0,324 

a 

-2n/5 

0.572 

0.543 

0.212 

0.218 

-0.490 

-0,468 

-0,361 

-0.373 

a 

-  n/2 

0.535 

0.498 

0.063 

0.064 

-0.604 

-0.570 

-0.417 

-0.433 

a 

-2n/3 

0./,88 

0.435 

-0.185 

-0.198 

-0.774 

•0.719 

-0.482 

-0.505 

a 

■  n 

0.424 

0.331 

-0.678 

-0.729 

-1.048 

-0.954 

-0.523 

-0.558 

We  remark  that  as  the  frequency  of  the  forced  oscillation  tends  to 
zero,  that  is  a  ♦  0,  Mybecomes  a  constant  and  being  a§*  the  damping 
factor  will  vanish.  In  two  dimensional  case  however,  being  O(loga)  the 
added  mass  tends  to  infinity  while  N^,  becomes  a  constant.  In  Figure  2,  as 
a  tends  to  zero  the  ordinate  becomes. 


ra 


-I  ♦  M  na  /.r2\  ^ 

3im  -jf-  N  --^ab  )  -  a 

a  ♦  o  ' 

where  A  -  nab,  and  ~  for  the  elliptic  disk.  Hence  at  the  origin 

a^ 


the  slopes  are, 


U3 


I 

1 

! 

I 

1 


■  9.8696 

for  b^l. 

lA  -  2.A67A 

for  ^1/A, 

1/8  -  1.2337 

for  b-1/8. 

respectively 

These  slopes  are  shown  In  Fifrure  2  by  the  stralf^ht  lines.  It  is  to  be 
observed  that  the  computed  results  deviate  very  rapidly  from  the  low 
frequency  approximation. 

For  the  roll  and  pitch,  the  results  comparable  to  those  presented 
in  Fi|o»re  3  to  FifJiure  6  cannot  be  found  in  the  existing  literatures. 

For  the  case  of  circular  disk  S»l,  realizing  that  roll  and  pitch  are  equi¬ 
valent  if  the  disk  is  turned  around  by  the  right  angle,  we  compare  tho  va¬ 
lues  of  fp(x,z)  and  f^(x,2)  for  roll  and  pitch  at  pivotal  points  of  the 
equal  radial  distance.  More  precisely  these  values  for  pitch  at  the  tip 
and  at  the  half  vadial  distance  on  Sf-axis  are  compared  with  those  for  roll 
at  the  tip  and  at  the  half  radial  distance  on  z-axis  in  Table  7  using  the 
results  of  tho  fine  lattice.  A  bad  agreement  was  found  at  a  ■  n  again  pro¬ 
bably  due  to  the  geometrically  unsuitable  lattice  in  use. 


Table  7 

Comparison  of  Real  and  Imaginary  Parts  of  Density  for  Roll  and  Pitch  at 
Pivotal  Points  of  Fx^ual  Radial  Distances  on  Circular  Disk, 

f_(5[.z)  f.(5c,z) 


At  Tip  on  3?,Tip  on  Z; 


(Roll) 

(Pitch) 

a 

■  n/6 

-I.IU 

-1.153 

a 

-  n/5 

-1.208 

-1.218 

a 

-  «//♦ 

-1.339 

-1.350 

a 

-  n/3 

-1.522 

-1.5<*5 

a 

-2r/5 

-0.905 

-0.977 

a 

■  n/2 

-0.126 

-0.1A6 

a 

-2n/3 

-0.265 

-0.268 

a 

-5«/6 

-0.5rtA 

-0.558 

a 

-  n 

-0.7/»8 

-1.026 

1/2  on  Sf, 

,1/2  on  Ti 

,Tip  on  X 

(Roll) 

(Pitch) 

(Roll) 

-0.669 

-0.657 

-0,020 

-0.739 

-0.722 

-O.OA6 

-0.886 

-0.859 

-0.138 

-1.193 

-1.150 

-0.671 

-0,795 

-0.812 

-1.259 

0.089 

0.060 

-o.euz 

0.259 

0.239 

-0.039 

0./.79 

O.IM* 

0.172 

0.561 

O.bOU 

-0.371 

Tip  on. 7, 

l/2  on  X, 

1/2  on  z 

(Pitch) 

(Roll) 

(Pitch) 

•0.020 

-0.012 

-0.011 

-0.0/U* 

-0.029 

-0.027 

-0.134 

-0.095 

-0.089 

-0.649 

-0.560 

-0.514 

-1.260 

-1.282 

-1.206 

-0.683 

-1.006 

-0.986 

-0.049 

-0.568 

-0.553 

0.209 

-0.396 

-0.378 

-0.215 

0.479 

0.276 

1 


Wf 

Vfe  further  observe  that  as  the  frequency  of  the  forced  oscillation 
tends  to  zero,  that  is  a  ^0,  I  and  I  become  constants,  while  H  and  H 
being  0(a*)  the  damping  factors  will  vanish. 

The  computation  work  is  performed  with  the  IBM  7090  data  processing 
system  at  the  Westinghouse  Electric  Corporation  in  East  Pittsburgh  Works, 
The  program  for  the  computation  is  coded  into  the  Fortran  language. 


Flrure  1 


1.0  r.o  3.0 

.1 
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1.0  r.o  :<.o 
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